Efficient classical-communication-assisted local simulation of n-qubit GHZ correlations 



Tracey E. Tessierjj Ivan H. Deutsch, and Carlton M. Caves 
Department of Physics and Astronomy, University of New Mexico, Albuquerque, New Mexico 87131 

(Dated: February 1, 2008) 

We present a local hidden-variable model supplemented by classical communication that repro- 
duces the quantum-mechanical predictions for measurements of all products of Pauli operators on 
an n-qubit GHZ state (or "cat state"). The simulation is efficient since the required amount of com- 
munication scales linearly with the number of qubits, even though there are Bell-type inequalities 
for these states for which the amount of violation grows exponentially with n. The structure of 
our model yields insight into the Gottesman-Knill theorem by demonstrating that, at least in this 
limited case, the correlations in the set of nonlocal hidden variables represented by the stabilizer 
generators are captured by an appropriate set of local hidden variables augmented by n — 2 bits of 
classical communication. 
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INTRODUCTION 

Bell's theorem Q codifies the observation that entan- 
gled quantum-mechanical systems exhibit stronger cor- 
relations than are achievable within any local hidden- 
variable (LHV) model. Beyond philosophical implica- 
tions, the ability to operate outside the constraints im- 
posed by local realism serves as a resource for many infor- 
mation processing tasks such as communication ^ and 
cryptography S]. 

The role of entanglement in quantum computation ^ 
is less clear, for the issue is not one of comparing quantum 
predictions to a local realistic description, but rather one 
of comparing a quantum computation to the efficiency 
of a realistic simulation. Nevertheless, various results in- 
dicate some connection between entanglement and com- 
putational power 1^, Q . Entanglement is a necessity if a 
pure-state quantum computer is to have scalable physi- 
cal resources 0- Moreover, systems with limited entan- 
glement can often be efficiently simulated classically 
Jozsa and Linden |^ showed that if the entanglement in a 
quantum computer extends only to some fixed number of 
qubits, independent of problem size, then the computa- 
tion can be simulated efficiently on a classical computer. 

Despite these results, global entanglement is by no 
means sufficient for achieving an exponential quantum 
advantage in computational efficiency [ijjj. The set of 
Clifford gates (Hadamard, Phase, and CNOT) acting on 
a collection of n qubits, each initialized to a fiducial state 
|0), can generate globally entangled states, yet accord- 
ing to the Gottesman-Knill (GK) theorem the out- 
comes of all measurements of products of Pauli operators 
can be simulated with 0{n^) resources on a classical 
computer. The GK theorem is an expression of proper- 
ties of the n-qubit Pauli group Vn \M i which consists of 
all products of Pauli operators multiplied by ±1 or ±i: 
the allowed (Clifford) gates preserve Vn, and the allowed 
measurements are the Hermitian operators in Vn- 

One approach to understanding the information pro- 



cessing capabilities of entangled states is to translate a 
quantum protocol involving entanglement into an equiv- 
alent protocol that utilizes only classical resources, e.g., 
the shared randomness of LHVs and ordinary classical 
communication. Toner and Bacon ,1^] showed that the 
quantum correlations arising from local projective mea- 
surements on a maximally entangled state of two qubits 
can be simulated exactly using a LHV model augmented 
by just a single bit of classical communication. Pironio 
took this analysis a step further, showing that the 
amount of violation of a Bell inequality imposes a lower 
bound on the average communication needed to repro- 
duce the quantum-mechanical correlations. 

In this article we analyze the classical resources re- 
quired to simulate measurements made on the n-qubit 
GHZ state ,14] (also called a "cat state"). We present 
a LHV model, augmented by classical communication, 
that simulates the quantum-mechanical predictions for 
measurements of arbitrary products of Pauli operators 
on this state. The simulation is efficient since the re- 
quired amount of communication scales linearly with n. 
These results are somewhat surprising in light of the ex- 
istence of Bell-type inequalities for n-qubit GHZ states 
where the amount of violation grows exponentially in the 
number of qubits [l5j |. 

Since the n-qubit GHZ state is generated by a circuit 
composed solely of Clifford gates, and since we consider 
only measurements of observables in Vn, our result yields 
an alternative perspective on the GK theorem. Whereas 
the GK simulation tracks the evolution of nonlocal hid- 
den variables that specify the generators of the stabi- 
lizer 0, 13, we simulate the circuit using local hidden 
variables that are supplemented by an efficient amount 
of classical communication to predict measurement out- 
comes. We conjecture that our result is general, i.e., that 
any GK circuit can be simulated with a LHV model plus 
an amount of communication that scales at most polyno- 
mially in the number of qubits. The existence of such an 
efficient classical model is currently under investigation. 
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FIG. 1: Circuit to generate the three-qubit GHZ state. 

Consider now the three-qubit GHZ state, = 
(|000) + |111)) /\/2, generated by the quantum circuit 
shown in Fig.^ In the language of the GK theorem, the 
evolution of the state is tracked by the evolution of the 
stabilizer generators. The Hadamard gate H transforms 
the Pauli operators X, F, Z according to HXH^ = Z, 
HYW = -Y, and HZW = X. Similarly, under the 
action of CNOT, we have 



C {XI) = XX 
C {YI) C'^ = YX , 
C {ZI) ^ZI , 



C {IX) = IX 
C {lY) = ZY , 
C{IZ)C^ = zz , 



(1) 



where the first qubit is the control, the second is 
the target, and / represents the identity operator. 
The stabilizer generators evolve through the circuit in 
Fig.mi^s {ZII, IZI, IIZ) ^^-^ civoTi. 

{XXI,ZZI,IIZ) ^^^^ 



> {XII, IZI, II Z) 

{XXX,ZZI,ZIZ). The 
full final stabilizer, consisting of all products of the gen- 
erators, includes -XYY, -YXY, -YYX, and XXX. 
This means that j-f/'s) is a -1-1 eigenstate of these four 
operators, which implies a deterministic violation of the 
assumptions of local realism . 

The GK description provides an efficient method for 
simulating the outcome of a measurement of any product 
of Pauh operators on the globally entangled state ji/'s), 
but it does so via the nonlocal stabilizer generators. We 
replace this nonlocal resource with a local description, 
augmented by classical communication, by constructing 
a LHV table where each row represents a qubit and each 
column represents a measurement. Locality is enforced 
by only allowing changes in rows corresponding to qubits 
that participate in an interaction. 

For the initial state |000), a measurement of Z on any 
qubit yields -1-1 with certainty, and a measurement of X 
or Y yields ±1 with equal probabilities. The first table in 
Fig. 121 gives corresponding LHVs for this state, with Rj 
denoting a classical random variable that returns ±1 with 
equal probability and j labeling the qubit to which the 
random variable refers. The table is read by choosing a 
measurement and multiplying the corresponding entries. 
The resulting product, with i discarded whenever it ap- 
pears, is the outcome predicted by the LHV model. The 
LHV table yields the correct quantum-mechanical pre- 
dictions for measurements of the 4^^ — 64 products of 
Pauli operators on the state |000). The use of i in the 
model, apparently just a curiosity, actually plays a cru- 
cial role. It simulates some of the conflicting predictions 
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FIG. 2: Evolution of the LHV model during creation of the 
three-qubit GHZ state. The rules for updating the LHV ta- 
bles are suggested by the equations for transforming Pauli 
operators. 



of commuting LHVs and anticommuting quantum oper- 
ators, which are the basis of Mermin's GHZ argument 
[l^. In addition, modehng the CNOT gates rehes on the 
X and Z entries being real and the Y entries being imag- 
inary and on the particular initial correlations between 
the X and Y values for each qubit. 

The first step in creating the three-qubit GHZ state 
is to apply the Hadamard gate to the first qubit. We 
extract rules for updating the LHV table from the trans- 
formations of the Pauli operators, which suggest that to 
simulate H, we should swap the X and Z entries, i.e., 
V° = Z'' and Z'' = X\ and flip the sign of the Y en- 
try, i.e., V° = —Y^, where b and a denote LHV values 
before and after a gate. Applying these rules to the first 
row leads to the second table in Fig. |21 which returns 
correct quantum-mechanical predictions for all measure- 
ments of Pauli products on the state (|0) + |1)) |00) /\/2. 
This is not surprising since the state remains a prod- 
uct state, and it is well known that a LHV model can 
be constructed for a single qubit 0. The usefulness of 
our model only becomes apparent when we apply it to 
entangled states. 

Applying the first CNOT gate in Fig. [T] yields the BeU 
entangled state (|00) |11)) |0) /\/2. To update the LHV 
table entries under a CNOT, we use the following rules 
for the control c and the target t: 



X- = X^ 



7b\rb 



Y- = Z'X 



b 7b 



(2) 



The update rules for U and CNOT keep the X and Z 
entries real and the Y entry imaginary, and the CNOT 
rule preserves the correlation XY Z = i that holds for 
each qubit after the operation of the Hadamard gate. 
Using the rules (0) in the first two rows gives the third 
table in Fig. |21to represent the Bell state. 

The LHV rules ^ must be consistent with the fif- 
teen transformations of nontrivial Pauli products under 
CNOT. For example, the transformation C (V/) = 
XX requires that Xl = X'^X^, which does follow from 
the rules (0). The CNOT rules (0) are derived from the 
six transformations , and because C = G\ these rules 
are consistent with five other transformations. Consis- 
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tency with the remaining four transformations, those be- 
ing C (XY) = YZ, C [XZ) Ct = -yy, and the in- 
verse transformations, requires that X^^^ = Y^Zi = 
Y^ZlZ\X\ and X\Z\ = -Y^Y^" = -Y^Z\X\Y^. These 
relations do not hold generally, but they are satisfied if 
the initial entries for both the control and target are cor- 
related according to XY Z — i, with X and Z real and Y 
imaginary. In all our applications of CNOT, these condi- 
tions hold. That they are not generally true is the chief 
obstacle to extending our results to arbitrary GK circuits 
and the entangled states they produce. 

The third table in Fig. [21 gives the correct quantum- 
mechanical predictions for all measurements of Pauli 
products on the BeU state (|00) |11)) |0) /\/2. What 
is new are the correlations between the first two qubits 
in each column. For example, the single-qubit measure- 
ments ZIl and IZI both return the random result 
the product of these outcomes always equals -1-1, the 
same as the outcome of a joint measurement of ZZl 
on the first two qubits. In this context, the j's in the 
correlated Y entries now lead to a problem: the single- 
qubit measurements YII and lYI both give the random 
result R1R2, with product -f 1, inconsistent with the out- 
come (iRiR2){iRiR2) = — 1 of a joint measurement of 
YYI [iTj. This problem persists throughout our anal- 
ysis, occurring for joint measurements involving K's on 
some qubits and having outcomes that are certain (i.e., 
measurements of stabilizer elements). This is the rea- 
son our LHV model must be supplemented by classical 
communication . 

At this point the problem is restricted to the joint 
measurements YYI and YYZ and the corresponding lo- 
cal measurements and thus can be corrected by flipping 
the sign of the outcome whenever a local measurement 
of Y is made on the first qubit; i.e., the model returns 
the random result —R1R2 for a measurement of YII. 
This sign flip fixes the required correlations and is ir- 
relevant to other joint measurements that involve Y on 
the first qubit, all of which have random results. Since 
the sign flip depends only on the measurement on the 
first qubit, it requires no communication between the 
qubits. Thus at this stage, with Bell-state entanglement, 
the LHV model gives correct quantum-mechanical pre- 
dictions for all observables in V3 and their correlations. 

We complete the simulation of the creation of the GHZ 
state by performing the CNOT between the first and 
third qubits, resulting in the last table in Fig.|21 This ta- 
ble yields correct quantum-mechanical predictions for all 
of the observables in 1^3, including those that form the 
basis of Mermin's GHZ argument flij, i.e., XXX = 1 
and XYY = YXY = YYX = -1. As promised, the 
imaginary Y entries make this agreement possible. 

Consider now the scheme for ensuring consistency with 
local measurement predictions for the three-qubit GHZ 
state. The only local measurements that yield inconsis- 
tent results are those associated with stabilizer elements 



that contain y's, i.e., the joint measurements XYY, 
YXY, and YYX. Let Alice, Bob, and Carol each pos- 
sess one of the qubits. If we put Alice in charge of en- 
suring compatibility, she should flip the sign of her out- 
come whenever she and/or Bob measures Y locally. This 
sign flip fixes the local correlations associated with XYY, 
YXY, and YYX and is irrelevant to other possible joint 
measurements that involve y's on the first two qubits, 
all of which have random outcomes. To implement this 
scheme, Bob must communicate to Alice one bit denoting 
whether or not he measured Y. For the three-qubit GHZ 
state, we thus have a LHV model, assisted by one bit of 
classical communication, that duplicates the quantum- 
mechanical predictions for all measurements in Vs and 
their correlations. 

The circuit that creates the general n-qubit GHZ 
state, lipn) = (|00. . .0) |11 . ..1))/V2, has the same 
topology as in Fig. ^ a Hadamard on the first qubit 
is followed by 71 — 1 CNOT gates, with the leading 
qubit as the control and the remaining qubits serv- 
ing successively as targets. The operator transforma- 
tions show that \ipn) is specified by the n stabilizer gen- 
erators (X^", . . . , 
The full stabilizer consists of the 2" observables in Vn 
that yield -1-1 with certainty. It contains Pauli products 
that have (i) only /'s and an even number of Z's and 
(ii) only AT's and an even number of Y^s, with an overall 
minus sign if the number of Y^s is not a multiple of 4. Of 
the remaining 2x4" observables in Vn , 2" are negatives 
of the stabilizer elements, thus yielding —1, while the rest 
return ±1 with equal probability ,4,]. 

Following the same procedure as in the three-qubit 
case, one finds that the LHV table representing the n- 
qubit GHZ state is given by 

X Y Z 

qubit 1 R2R3 ■■■Rn iRiR2 ■ ■ ■ Rn Ri 
qubit 2 i?2 iRiR2 Ri 

qubit 3 R3 iRiRs Ri ■ (3) 

qubit n Rn iRiRn Ri 

That this table gives the correct quantum-mechanical 
predictions for all measurements of Pauli products fol- 
lows from the consistency of our LHV update rules, but 
it is nevertheless useful to check this directly. Suppose a 
Pauli product contains no A"'s or Y's, but consists solely 
of /'s and Z's. Then it is clear from the table in Eq. Q 
that the outcome is certain if and only if the number of 
Z's in the product is even. Suppose now that the product 
has an A" or a y in the first position. Then it is apparent 
that to avoid a random variable in the overall product, 
all the other elements in the product must be X's or Y's 
and the number of Y's must be even; the outcome is +1 
if the number of K's is a multiple of 4 and —1 otherwise. 
Finally, suppose the Pauli product has an A or a y in a 
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position other than the first. Then the only way to avoid 
a random variable in the overall product is to have an 
X or a y in the first position, and we proceed as before. 
This argument shows that the LHV table for the n-qubit 
GHZ state gives correct quantum-mechanical predictions 
for measurements of all Pauli products. 

It remains to ensure that the products of the LHV pre- 
dictions for local measurements are consistent with the 
corresponding joint measurement results. As before, the 
source of the inconsistency is the i in the Y table entries, 
the very thing that allows us to get all the Pauli products 
correct. Stationing Alice at the first qubit and putting 
her in charge of ensuring consistency, we see that what 
she needs to know is the number of i's in the product for 
the corresponding joint measurement. In particular, let- 
ting qj = i if y is measured on the jth qubit and qj — 1 
otherwise, Alice can ensure consistency by changing the 
sign of her local outcome if the product p„ — qi ■ ■ ■ qn 
is —1 or — i and leaving her local outcome unchanged if 
p„ is or j. This scheme requires n — 1 bits of com- 
munication as each of the other parties communicates to 
Alice whether or not he measured Y, but we can do a bit 
better. Alice's action is only important when p„ is -1-1 
or —1; when p„ is i or — i, the sign flip or lack thereof 
is irrelevant because the joint measurement outcome is 
random. As a result, Alice can get by with the truncated 
product Pn-i = 51 • • • qii-i- she flips the sign of her local 
outcome if p„_i is i or — 1 and leaves the local outcome 
unchanged if is — i or 1. The scheme works be- 

cause whether qn is 1 or i, Alice flips when p„ = — 1 and 
doesn't flip when p„ = -1-1, as required. This improved 
scheme requires n — 2 bits of classical communication; it 
generalizes our previous results for the Bell state and the 
three-qubit GHZ state. 

The consistency scheme generalizes trivially to the case 
of measurements made on I disjoint sets of qubits. For 
each set k chosen from the I sets, the table yields a mea- 
surement product that is the predicted outcome multi- 
plied by qk = i or qk = 1. Letting Alice be in charge of the 
flrst set, all but the last of the other sets communicates 
qk to Alice, who computes the product qi - ■ ■ qi-i and 
decides whether to flip her set's outcome just as before. 
Consistency with the corresponding joint measurement is 
thus ensured at the price of / — 2 bits of communication. 

Using local hidden variables and an efflcient amount of 
classical communication, we have shown that it is possi- 
ble to simulate the correlations that arise when measur- 
ing arbitrary products of Pauli operators on an n-qubit 
GHZ state. Though the n-qubit GHZ state is highly 
entangled, the probability distributions for the allowed 
measurements of Pauli products are essentially trivial, 
being either certainty or binary randomness. This prop- 
erty is shared by all states produced by GK circuits, lead- 
ing us to conjecture that our results can be extended to 
measurements of Pauli products on any state produced 



by a GK circuit. In contrast, allowing just one addi- 
tional nontrivial measurement, say of {X -\- Z) /v^, leads 
to correlations for which our simple simulation will no 
longer work. We anticipate that under this more gen- 
eral measurement scheme, the amount of classical com- 
munication required to make a LHV model work grows 
exponentially in the number of qubits. 

Our model provides weak evidence that the power of 
quantum computation arises not directly from entan- 
glement, but rather from the lack of an efflcient, lo- 
cal realistic description assisted by an efflcient amount 
of nonlocal, but classical communication. An efficient 
communication-assisted LHV model for all GK circuits 
would provide powerful additional evidence for this idea. 
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